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Abstract 

String theory suggests the simultaneous presence of many ultralight axions, possibly 
populating each decade of mass down to the Hubble scale 10~ 33 eV. Conversely the presence 
of such a plenitude of axions (an "axiverse" ) would be evidence for string theory, since it arises 
due to the topological complexity of the extra-dimensional manifold and is ad hoc in a theory 
with just the four familiar dimensions. We investigate how several upcoming astrophysical 
experiments will be observationally exploring the possible existence of such axions over a vast 
mass range from 10 _33 eV to 10~ 10 eV. Axions with masses between 10 _33 eV to 10~ 28 eV can 
cause a rotation of the CMB polarization that is constant throughout the sky. The predicted 
rotation angle is independent of the scale of inflation and the axion decay constant, and is of 
order a ~ 1 /137 -within reach of the just launched Planck satellite. Axions in the mass range 
lCT 28 eV to lCT 18 eV give rise to multiple steps in the matter power spectrum, providing us 
with a snapshot of the axiverse that will be probed by galaxy surveys-such as BOSS, and 21 
cm line tomography. Axions in the mass range 10 _22 eV to 10~ 10 eV can affect the dynamics 
and gravitational wave emission of rapidly rotating astrophysical black holes through the 
Penrose superradiance process. When the axion Compton wavelength is of order of the black 
hole size, the axions develop "superradiant" atomic bound states around the black hole 
"nucleus". Their occupation number grows exponentially by extracting rotational energy 
and angular momentum from the ergosphere, culminating in a rotating Bose-Einstein axion 
condensate emitting gravitational waves. For black holes lighter than ~ 10 7 solar masses 
accretion cannot replenish the spin of the black hole, creating mass gaps in the spectrum of 
rapidly rotating black holes that diagnose the presence of destabilizing axions. In particular, 
the highly rotating black hole in the X-ray binary LMC X-l implies an upper limit on the 
decay constant of the QCD axion f a < 2 x 10 17 GeV, much below the Planck mass. This 
reach can be improved down to the grand unification scale f a < 2 x 10 16 GeV, by observing 
smaller stellar mass black holes. 



The Principle of Plenitude: "This best of all possible worlds 
will contain all possibilities, with our finite experience of 
eternity giving no reason to dispute nature 's perfection. " 



Gottfried Leibniz, "Theodicee" . 
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1 A Plenitude of Axions 
1.1 The QCD Axion 

The Standard Model QCD action can be appended by the CP-violating topological interaction p] 

Se= 3^J d A xe^Ti G, u G Xp . (1) 

This term is a total derivative and does not contribute to the classical field equations. However, it 
affects physics at the quantum level due to the existence of topologically non-trivial field configu- 
rations. Under the shift 9 — ► 9 + 2n the action changes by 2n, indicating that 9 is a periodic 
parameter, with a period equal to In . In the presence of fermions, as a consequence of the chiral 
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anomaly, the actual physical parameter is not 9 itself, but 9 — the sum of 9 and the overall phase 
of the quark mass matrix 



Stringent bounds on the neutron electric dipole moment, d n < 2.9 ■ lCr 26 ecm [2], imply that 9, if 
non-zero, should be tiny 9 < 1CT 10 . This is the strong CP problem — one of the most tantalizing 
hints for physics beyond the SM. It shares some important features with other more severe fine- 
tuning problems of the SM — the cosmological constant and the weak hierarchy problems. Just 
as its more famous cousins, the strong CP problem requires an extreme fine-tuning of apparently 
unrelated quantities, 9 and argdetm ? , to produce a tiny 9 as required by observations. Moreover 
because of the observed CP-violation in the quark sector no symmetry of the SM is restored 
as 9 — > 0. The important difference with the cosmological constant and hierarchy problems is 
that the smallness of 9 does not appear to be a necessary prerequisite for the existence of life. 
Consequently, unlike these other mysteries there is no ambiguity on whether the strong CP problem 
has an anthropic or a dynamical solution. The smallness of 9 is a clear call for a new dynamical 
mechanism. 

An attractive solution to the strong CP problem is to promote 9 to a dynamical field — the 
QCD axion a. At the classical level the axion action is assumed to be invariant under shifts 
a — ► a + const. In other words, at the classical level the axion is a Nambu-Goldstone boson of 
a spontaneously broken global symmetry (Peccei-Quinn (PQ) symmetry) |3HH|3]. This symme- 
try is preserved in the quantum theory at the perturbative level. However, the PQ symmetry is 
anomalous in the presence of topologically non-trivial QCD field configurations, explicitly break- 
ing the shift symmetry and generating a periodic potential for the axion. The resulting axion 
vacuum expectation value (vev) then automatically adjusts itself to cancel 9, solving the strong 
CP problem. The physical properties of the QCD axion are to large extent determined by the 
scale f a of the PQ symmetry breaking, similar to how the low energy pion interactions are fixed 
by the pion decay constant /„.. In particular the axion mass m a is given by 



Apart from the coupling to gluons ([I]) the QCD axion may have similar couplings to other gauge 
bosons, most notably photons, and derivative couplings to fermions. The precise numerical values 
of these couplings are model dependent, but their overall scale is determined by the axion decay 
constant f a . Laboratory searches and astrophysical constraints exclude values of f a below ~ 10 9 
GeV and, for sufficiently large couplings to photons and, assuming dark matter abundance for the 
QCD axion, a small region around 10 11 GeV [6]. 

It would be especially interesting to find a QCD axion with f a ^> 10 12 GeV. Indeed, in this 
case for typical initial values of the axion field the relic abundance of axions would be higher than 
the critical density leading to the overclosure of the Universe El E] . This does not imply that 
high faS are excluded — if the observed (quasi)homogeneous Universe were just a small patch of 
the inhomogeneous Multiverse, then the axion initial conditions may be different from place to 
place and life can only develop in rare regions with an atypically small axion density [TU]. This 
anthropic argument is of a particularly mild form — the post-inflationary probability distribution 



9 = 9 + arg det m q . 



(2) 




(3) 
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of the axion field in the Multiverse is calculable since the axion properties are entirely encoded in 
a low-energy effective quantum field theory and the axion is defined in the finite range (0, 2irf a ) 
Consequently, finding a high f a QCD axion would not only explain the smallness of 9, but 
also provide strong support to the idea of eternal inflation, which is a natural way to create the 
Multiverse, and possibly the idea of the string landscape more generally. 

1.2 String Theory Axions 

We see that the axion provides a potentially testable dynamical solution to the strong CP problem 
and its discovery may have even more profound implications. However, at the effective field theory 
level it is hard to judge how natural it is to have such a "fake" global PQ symmetry which is 
explicitly broken exclusively by QCD. Note, that in order not to spoil the solution to the strong 
CP problem all other sources of explicit PQ symmetry breaking should be at least 10 orders of 
magnitude down with respect to the potential generated by the QCD anomaly, and one may 
be especially cautious about the viability of such a proposal, given the common lore that global 
symmetries always get broken by quantum gravitational effects [121 H31 EH] ■ This makes it natural 
to inquire whether axions arise naturally in the most developed quantum theory of gravity — string 
theory. 

Pseudoscalar fields with axion-like properties generically arise in string theory compactifica- 
tions as Kaluza-Klein (KK) zero modes of antisymmetric tensor fields [T5l [16] . Examples of such 
fields include, for instance, the Neveu-Schwarz 2-form B2, present in all string theories, or the 
Ramond-Ramond forms 60,2,4 of type IIB theory, or C13 of type IIA theory. An interesting prop- 
erty of higher-order antisymmetric tensor fields as opposed to scalar fields (0-forms), is that upon 
compactification they typically give rise to a large number of KK zero modes, determined by the 
topology of the compact manifold M.§. For instance, the number of massless scalar fields resulting 
from a single two-form Bmn or Cmn is equal to the number of (homologically non-equivalent) 
closed two-cycles Cj in A4q. Namely, by taking both of the indices of the two-form along a cycle one 
gets a massless (pseudo) scalar four- dimensional field corresponding to that cycle. More formally, 
the KK expansion for the two-form contains 



where x are non-compact coordinates, y are compact coordinates, and uji form the basis of closed 
non-exact two-forms (cohomologies) dual to the cycles C iy i.e., satisfying 



A similar argument applies to higher-form fields, for instance the number of (pseudo)scalar zero 
modes corresponding to C4 is equal to the number of homologically non-equivalent four-cycles. Of 
fundamental importance to what follows is the fact that in the vast majority of compactifications 
the number of cycles is quite large, of order several hundred, even up to 10 5 [T7j. The basic reason 
for this is simple combinatorics — it is natural to expect many non-equivalent embeddings of a 
closed two-surface into a reasonably complicated six-dimensional manifold. For instance, even 




(4) 




(5) 
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the simplest Calabi-Yau manifold — the six-torus — has (6 x 5)/2 = 15 different two-cycles and the 
same number of four-cycles. Of course, this example is too simple to give rise to a realistic gauge 
and matter field spectrum. Typically, even the simplest string compactifications incorporating 
the Standard Model gauge group and chiral fermion content give rise to more than a hundred 
two- cycles. 

The four-dimensional scalar fields resulting from this KK reduction are not only massless, but 
have zero potential as a consequence of the higher-dimensional gauge invariance of the antisym- 
metric tensor field action. This invariance guarantees that no potential is generated at any order 
of perturbation theory. However, antisymmetric tensor fields have Chern-Simons couplings, which 
are crucial for the Green-Schwarz mechanism for anomaly cancelation. Upon KK reduction these 
terms can result in the axionic couplings <Q to gauge fields [H] . In type IIB theory this happens, 
for instance, for a C<i axion if there is a D5 brane wrapped over the corresponding two-cycle, or 
for a C4 axion if there is a D7 brane wrapped over the corresponding four-cycle, both of which 
are natural occurrences in the string landscape. So we see that string theory has a potential of 
providing many particles with the qualitative properties of the QCD axion. 

Nevertheless, the above arguments do not suffice to say that the QCD axion is predicted by 
string theory. First, string axions can be at tree level completely removed from the spectrum of 
light fields by the presence of fluxes, branes and/or orientifold planes [T71 [T5]. For instance, the 
DBI action of a D5 brane depends on the NS two-form B MN , so that wrapping such a brane around 
a two-cycle makes the corresponding axion heavy, with mass of order the string scale. Similarly, 
the RR C*2 axion can be lifted by wrapping a NS5 brane over the corresponding two-cycle. From 
now on we will focus exclusively on the axions that escape such tree-level liftings. 

Secondly, even if an axion does not become heavy due to these tree-level effects its potential 
always acquires non-perturbative contributions from one or more of a variety of sources. These 
include, worldsheet instantons [20], euclidean D-branes wrapping the cycle [21], gravitational 
instantons [2], and gauge theory instantons [1] if the axion couples to a non-abelian gauge group. 
In many cases these corrections can be large enough to ruin the solution of the strong CP problem. 

We see, that one cannot really make a case that the QCD axion is predicted by string theory, 
but rather the requirement for one of the string axions to be responsible for a solution of the 
strong CP problem puts a restriction on string theory model building. For instance, it disfavors the 
possibility that all string moduli are stabilized by a supersymmetry(SUSY)-preserving mechanism, 
as happens, for instance, in the KKLT scenario [22]. Indeed, if that were the case axions would 
receive large SUSY preserving masses together with their superpartners — saxions. On the other 
hand if not all of the moduli are stabilized in a supersymmetric way, saxions (and axinos) receive 
masses after SUSY breaking, while axions, being (pseudo)Nambu-Goldstone bosons, are protected 
from these contributions. 

Note that SUSY breaking masses for moduli typically come out in the 1 GeV -hi TeV range. 
This gives rise to the infamous cosmological moduli problem — light scalar moduli may overdose 
the Universe and/or their late decays may destroy successful BBN predictions and introduce 
non-thermal distortions into the CMB spectrum. The most straightforward way to get around 
this problem is to assume that the expansion rate during inflation was relatively low, such that 
the moduli are heavy during inflation and not produced, Hi n fi < 0.1 GeV. The corresponding 
inflationary energy scale is still quite high Ei n fi ~ 10 8 GeV, so that if reheating is efficient the 



4 



temperature is high enough for successful baryogenesis. Such a cosmological scenario will be 
assumed throughout this papeiQ 

As we will argue now, even though one cannot really predict the existence of the QCD axion in 
string theory, the assumption that the QCD axion does exist and is responsible for the smallness of 
9, combined with the facts about string theory axions reviewed above, not only puts restrictions 
on string theory model building, but also strongly suggests a rather predictive scenario with a 
distinctive set of signatures for upcoming cosmological and astronomical observations. 

The point is, as we already stressed, there is no reason for the solution to the strong CP 
problem to have an anthropic origin. Consequently, neither the very existence of the QCD axion, 
nor the extreme smallness of all non-QCD contributions to its potential should be a result of 
a fine-tuning. Instead, it has to be a natural dynamical consequence of the properties of the 
compactification manifold giving rise to the string theory vacuum where we live. But given that 
string theory compactifications have a potential of producing hundreds of axions it would be really 
surprising in such a situation if only the QCD axion were then light. Note also that of all axions 
only one linear combination gets a mass from the explicit breaking due to QCD. Consequently, we 
come to the conclusion that if one takes seriously the QCD axion as a solution to the strong CP 
problem then in string theory one expects to find many light axions. As we will see, these axions, 
if they exist, can be observed in a number of different ways. 

The two principal parameters characterizing a general string theory axion are its mass m and 
decay constant f a . Unlike for the QCD axion these two parameters are not related by What 
values for these parameters can one expect? First, it's convenient to parametrize the effective 
four-dimensional axion Lagrangian as 

C = f l (da) 2 - A 4 f/(a) , (6) 

where U(8) is a periodic function with a period equal to 2ir. The overall scale of the potential A 4 
is related to f a and m through 

m = A 2 /f a . (7) 

Since it is non-perturbative effects that provide the potential in cases where the axion escapes tree 
level lifting, we write the scale in the following form, 

A 4 = A" 5 , (8) 

where \i is a UV energy scale. In general more than one type of non-perturabtive effect contributes 
to the potential so K A U (a) should be thought of 

i 

where each of the dynamical scales A 4 has the form If the axion potential arises from the 
superpotential generated by string instantons, then the UV scale /i is the geometric mean of 

1 In principle this can be relaxed if there is no low energy SUSY. Note, that still there is a quite restrictive bound 
on the inflation scale from isocurvature perturbations [23 , 24 , if the PQ symmetry is broken during inflation, which 
is the natural option for the string theory axions. 
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the string/Planck scale (which sets the scale of instanton physics in this case) and of the SUSY 
breaking scale Fslsy, 

// 4 ~ F susy M 2 Pl . (10) 
In string theory constructions one finds, in many cases, that f a and S are related by [25j dEJ [26] 

f Mpi rm 

This relation has a simple geometrical origin. Namely, the four dimensional Planck mass is 
determined by the typical size L of the compactification manifold as 

M 2 Pl ~ 97 2 LX 8 , (12) 

where l s and g s are the string length and the string coupling. Similarly, the axion decay constant 
is also determined by these parameters and by the area A of the cycle. For instance, for an axion 
coming from the RR two-form in IIB theory integrated on a 2-cycle of area A ~ L 2 one has 

f a ~ L%*A- 2 ~ IV , (13) 

where the two factors of A~ x came from the powers of the inverse metric in the kinetic term of the 
two-form before KK reduction. The action of the string instanton generating the axion potential 
is given by the product of the tension of the euclidian D brane that wraps the cycle and the area 
of the cycle, 

S ~ g- l K 2 L 2 , (14) 
thus giving the relation ([TTj)J^| The value of f a can be significantly lower than the esimate (11) 



if a compactification manifold is significantly anisotropic, or, particularly, if a large amount of 
warping is present. However, in explicit examples f a never comes out parametrically higher than 



(11) and it was argued that the inequality 



fa < (15) 

follows from very basic principles of black hole physics ("gravity as the weakest force" conjecture) 
Throughout most of the discussion in this paper we assume that the compactification man- 



ifold is not too anisotropic and the amount of warping is limited, so that the relation (11) is a 
good guide for the possible values of f a . 

In the case of the QCD axion not only QCD non-perturbative effects but also string instantons 
give contributions to the potential. To solve the strong CP problem and not lead to a too large 9, 
string instanton contributions to the potential must be subdominant by a factor of 10 10 compared 
to QCD, implying, for intermediate scale supersymmetry breaking, a string instanton action 

S > 200 . (16) 



2 Similarly, if the axion potential is generated by the strong dynamics of a gauge sector coming from a stack of 
D5 branes wrapping the two-cycle, the gauge instanton contribution goes as e~ s with S ~ 2tt /a where the inverse 



gauge coupling at the compactification scale is given by 1/a ~ L 2 /g s l 2 s so that one again arrives at (111]). Clearly, 
all these relations hold up to order one numerical factors. 
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This discussion then suggests the following scenario for the distribution of f a and m for different 
axions. The values of f a are inversely proportional to the area of the corresponding cycle, so they 
do not change much from one axion to another. Given that the compactification is such that 
S > 200 for string contributions to the QCD axion, and no special fine tuning is allowed, all axion 
decay constants in this scenario are likely to be close to the GUT scale Mqut — 2 x 10 16 GeV. 
On the other hand, axion masses are exponentially sensitive to the area of the cycles, so that 
we expect their values to be homogeneously distributed on a log scale. Given that, as argued 
above, one can expect several hundred different cycles this suggests that there may be several 
string axions per decade of energy. It has also been argued recently that the mixing of axions 
with vacuum energy cancelling 4-forms in the Bousso-Polchinski landscape can yield axion mass 
scannings [28]. In Section 2 we show that there are a variety of observational windows on such 
axions. 

1.3 Wilsonian Scanning of the Cosmological Constant 

Interestingly, the same number — of order several hundred cycles — is also suggested by the require- 
ment that there are enough different fluxes to implement the Bousso-Polchinski mechanism for 
scanning the cosmological constant [29, 30] . The famous 10 500 estimate for the number of different 
string vacua comes out exactly as a result of having 500 different fluxes, corresponding to different 
cycles. At this point one may wonder whether the axion potential itself may be responsible for 
the scanning of the cosmological constant in the presence of several hundred axions. 

Indeed, as explained in Ref. [31], in the presence of N decoupled scalar fields with two or more 
non-degenerate minima each with energies of order Mq UT there are of order 2^ minima that scan 
the total vacuum energy down to ~ 2~ N Mq UT . In fact, there is no reason why the minima of the 
individual fields should all be at the same scale Mq UT . For instance, the opposite extreme is to 
have at each energy scale a field with several minima, which cancels the vacuum energy at this 
particular scale. An intuitive way to think about this kind of scanning is to write the vacuum 
energy in the binary system. If each field has two minima, the choice of the minima amounts to 
putting either or 1 at the corresponding digital place. We will refer to this model of scanning 
as Wilsonian scanning. Of course, in a generic case, instead of one of these two extreme cases 
of either uniformly distributed scalars or all scalars at the GUT scale, there could be densely 
populated mass ranges with gaps in between them. The general condition for successful scanning 
is that the number of scalars above any energy E is larger than ~ log(Mpi/E ). 

In order that axions be responsible for such scanning, it is necessary that the axion potential has 
non-degenerate minima. If the axion potential is generated by a string instanton it is dominated by 
a single contribution (because the ra-instanton action is nS) and the periodic function U (9) in ^ 
is simply U(8) = 1 — cos 9 and there is only a single minimum. On the other hand, if the potential 
is generated by strong IR gauge theory dynamics there is no reason for a strong suppression of 
the higher instanton contributions, so U(6) is a generic periodic function and it is natural for it to 
have several non-degenerate minima. A simple toy example of such a situation is when the one- 
and two-instanton contributions enter with comparable coefficients, U(9) ~ 2 — (cos 6* + cos 26). 

Consequently, if light axions are (partially) responsible for the scanning of the cosmological 
constant at low energies they should be accompanied by a large number of strongly coupled QCD- 
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like hidden sectors with low confinement scales. This both opens up interesting phenomenology 
associated to the presence of this "dark world" and raises the question of how it managed to escape 
being observed so far. We will touch on some of the issues involved in the concluding Section [3} 
For now we focus upon the observational signatures of the light axions that we have argued are 
generic to string theory once the strong CP problem is solved. 



2 Cohomologies from Cosmology 
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Figure 1: Map of the Axiverse: The signatures of axions as a function of their mass, assuming 
f a ~ Mqut and Hi n f ~ 10 8 eV. We also show the regions for which the axion initial angles are 
anthropically constrained not to over-close the Universe, and axions diluted away by inflation. 
For the same value of /„ we give the QCD axion mass. The beginning of the anthropic mass 
region (2 x 10~ 20 eV) as well as that of the region probed by density perturbations (4 x 1CT 28 
eV) are blurred as they depend on the details of the axion cosmological evolution (see Section 
3). 3 x 10~ 18 eV is the ultimate reach of density perturbation measurements with 21 cm line 
observations. The lower reach from black hole super-radiance is also blurred as it depends on 



the details of the axion instability evolution (see Section 2.5). The region marked as "Decays 



outlines very roughly the mass range within which we expect bounds or signatures from axions 
decaying to photons, if they couple to E-B. We will discuss axion decays in detail in a companion 
paper. 



2.1 Discovering the String Axiverse 

We now turn to the observational consequences of axions lighter than or around the QCD axion 
mass. For simplicity, we keep f a fixed at Mqut and H in fi ~ 0.1 GeV. The initial displacement of 
axions heavier than ~ 10~ 20 eV has to be tuned in order for them not to overdose the universe and 
axions heavier than 0.1 GeV have been diluted away by inflation. The observational consequences 
of the string axiverse are outlined in Figure [TJ 



We concentrate on three main windows to the axiverse. First, as discussed in Section 2J2 
axions of masses between 10~ 33 eV and 4 x 10 -28 eV, if they couple to E ■ B, cause a rotation in 
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the polarization of the Cosmic Microwave Background (CMB) by an angle of order 10~ 3 , which is 
close to the current bounds (~ 10~ 2 ) [32], [33] . Future experiments, such as Planck and CMBPol 
will be able to probe values of the rotation in the CMB polarization down to 10~ 5 [34J. 

Second, as discussed in Section 2J3 axions with masses higher than 10~ 28 eV can be a significant 
component of the dark matter (DM) and suppress power in small scale density perturbations 
(< 1 Mpc). This is because the quantum pressure scale originating from the uncertainty principle 
and below which gravitational collapse is not possible, is proportional to 1/ \fm and thus for these 
light axions is a cosmologically observable scale. Since the axiverse should contain a plethora of 
string axions with masses homogeneously distributed on a log scale, the existence of multiple steps 
in the small scale perturbation spectrum is a natural expectation. The amount of the suppression — 
the step height — is proportional to the fraction of the DM constituted by the particular axion. 
Such steps may be detectable with the BOSS [35] and 21 cm line observations [36]. In particular, 
the 21 cm line tomography will be sensitive to masses up to 3 x 1CT 18 eV that are well inside the 
anthropic regime. 

Finally, axions of masses between 10~ 22 and 10~ 10 eV can affect the dynamics of rotating black 
holes due to the effect of superradiance. When a black hole rotates, a boson with a Compton 
wavelength comparable to the black hole size creates an exponentially growing bound state with 
the black hole. This gravitational atom can be de-excited through graviton emission that carries 
away the black hole's angular momentum. For black hole masses larger than ~ 10 7 M , or axion 
masses smaller than 1CT 18 eV, accretion may still be efficient enough to support the maximal 
rotation and sustain a "Carnot cycle" that turns the black hole into a gravity wave pulsar with 
possibly detectable signal at future gravity wave experiments. For lighter black holes (heavier axion 
masses) this effect leads to a spindown of the black hole, resulting in gaps in the mass spectrum 
of rapidly rotating black holes. With the quality of data constantly improving, measurements of 
the spin of stellar mass (~ 2 4- 10 M Q ) black holes will be able to probe also the QCD axion 
parameter space with f a > 10 16 GeV, well inside the region where the QCD axion relic abundance 
is anthropically constrained. These effects are discussed in Sections 2.4 and 2.5, and in the 
Appendix. 

For axions in the range ~ 10~ 9 to ~ 10 s eV, and assuming the axions have couplings to E ■ B, 
decays to photons can potentially lead to signatures. A companion paper will discuss such decays, 
as well as the physics induced by warping the axion decay constant to scales lower than Mqut, 
and the many dark sectors implied by Wilsonian scanning and/or highly warped throats. 



2.2 Rotation of the CMB Polarization 

Axions much lighter than the QCD axion, when they have an E ■ B coupling to electromagnetism 
(EM), change the polarization of the CMB photons if they start oscillating anytime between 
recombination and today. These axions cannot couple to QCD, as 4d gauge coupling unification 
implies, otherwise they would get large contributions to their masses. A coupling to QCD, however, 
can be easily avoided in the framework of orbifold GUTs [37] [3S] EH HU] • An example of such 
a theory is a scenario with one extra dimension where SU(5) is preserved in the bulk and the 
breaking down to SU(3) C x SU(2)l x U(1)y occurs on the boundary of the extra dimension. The 
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SM gauge couplings are given by 



\ = ^ + T2> ^thc =1,2,3 (17) 

9 C 95 n c 

where V is the extra-dimensional volume, g 5 is the 5d SU(5) coupling and h c are the gauge 
couplings of the SU (3) x SU(2) x U(l) brane kinetic terms that are allowed by the brane-localized 
breaking of SU(5). When the volume of the extra dimensions is parametrically large (equivalently 
when the effective 4d coupling of the SU(5) is much smaller than the brane-localized gauge 
couplings), there is apparent SU(5) unification for the SM gauge couplings, 

withe =1,2, 3, (18) 

9l 91 

with corrections that are parametrically of the same size as traditional GUT-scale threshold cor- 
rections. In these scenarios, since gauge coupling unification is not true everywhere in the extra 
dimensions, axions that have brane localized couplings, naturally couple to SU(2) L or U(l)y 
without coupling to QCD. 

Only a few "local" axions are able to couple to electromagnetism with full strength in this 
way. Most axions resulting from antisymmetric forms are localized on cycles far away from the 
position of the SM in the full compactification. On the other hand, axions from cycles intersecting 
with ours have a kinetic mixing with our axion and could be more weakly coupled to SU(2)l or 
U(1)y- Specifically, defining 

7# = / u>i A* Uj , (19) 
Jm 

where tOi is the basis of the closed two-forms dual to the cycles C\ as in Q, the kinetic terms of 
the four-dimensional axion fields are of the form 



d'xji^bidfjbj . (20) 



The axions also receive a variety of contributions to their potentials leading to spectrum of axion 
masses m^. The end result after diagonalization is that the axions from intersecting cycles also 
acquire a coupling to electromagnetism, which is suppressed, however, by the mixing angle 0y ~ 
mj/m 2 j between the axions. This can be a significant suppression for widely separated dynamical 
scales Aj. 

For axions that couple only to EM, interactions are summarized by the following Lagrangian 
£ = + \(d^) 2 - KU(Hfa) + ^e^F,„F Xp . (21) 

Here a dimensionless constant C is of order one for an axion that coupled to the photon directly. 
For instance, in the four-dimensional SU(5) GUT it is equal to C — 4/3. When the photon wave- 
length is short compared to a typical length scale of variation in the axion field, the combinations 

D = (e+ \ ^y$B^ and H = (b — ^^j-(f)E^j satisfy free wave equations [H] . Consequently, a 
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region of space with an inhomogeneous axion field becomes optically active [42| |4"1T |4*3| HI] — a 
linearly polarized freely propagating photon experiences a rotation of the polarization plane by 
an angle A/3 equal to 

W = ^r [dr<P, (22) 

^Ja J 

where r is the time along the photon trajectory. Due to frequent Compton scatterings off electrons 



CMB photons are not polarized before recombination, so the integration in (22) goes from the 
time of recombination, r rec , till today, roj^] As a result one obtains for the rotation angle 

A/3 = ^ (0(r o ) - 0(r rec )) . (23) 

Note that this result is only valid for sufficiently small axion masses, m < St~^ c , where St rec ~ lOkpc 
is the duration of recombination. For larger masses rapid oscillations of the axion field on the 
timescale of recombination lead to suppression of the rotation angle by a factor e _m<5 ' rec . The 
rotation angle is maximum for axion masses smaller than the expansion rate at recombination 
Hrec ~ (7 x 10 5 yr) 1 and bigger than the current Hubble parameter Hq ~ (10 10 yr) , while it 
is still substantial for masses up to ~ 10H rec . In this regime the axion field has negligible value 
today, 0(to) ~ 0. On the other hand at recombination axion oscillations had not yet started, so 
that the axion field took its primordial value set during inflation, |0(r rec )| ~ irf a /\/3. As a result, 
the rotation angle from axions of mass between 10~ 33 eV and 4 x 10~ 28 eV becomes, 

AP ~ 2^1 W 10_3 ' (24) 

If iV axions are present in this mass range, the rotation angle gets enhanced by a factor of 
y/N. The remarkable feature of this result is that the rotation angle (24) is independent of both 
the PQ scale f a , and the scale of inflation. As previously mentioned, we assume that H in f ~ 0.1 
GeV, so isocurvature fluctuations from the axion field are extremely suppressed. Consequently, 
the rotation angle A/3 is constant throughout the sky. As a result of this rotation a part of the 
E-mode polarization gets converted into the B-mode. Therefore, even though, as a consequence of 
a low inflationary scale, we do not expect the gravitational wave B-mode signal in the axiverse, B- 
mode can nevertheless be generated through this effect. The two sources of the B-mode are easily 
distinguishable, because axions generate also EB and TB cross-corellations which are forbidden 
by parity in the standard case. 

Even though current limits [32J, [33] on the rotation angle are at the level A/3 < 2° = 3.5 • 10~ 2 , 
i.e. above the level of the expected signal for one axion, Planck will improve this limit by an order 



of magnitude down to 0.1° and will be sensitive to (24). Moreover it is expected that CMBPol 
will be sensitive to a signal as low as 0.005° [31] . Consequently, in the future we will be sensitive 
to axions even somewhat outside the optimal mass range, or to smaller axion couplings to E • B, 



when the signal is suppressed as compared to (24). 



3 Photons also experience Compton scattering during reionization. To incorporate this effect one has to solve 
the appropriate kinetic equations rather than just the Maxwell's equations with an axionic source. As soon as the 
axion mass is much larger than the Hubble parameter of the Universe at reionization this effect is unimportant and 



the integral in (22 1 is saturated at early times. Also it affects the spectrum mainly at small Vs, where the cosmic 
variance is large. 
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2.3 Steps in the Power Spectrum 

A purely gravitational signature of light axions is the presence of step-like features in the matter 
power spectrum. The suppression of the CDM power spectrum at small scales in the presence 
of ultra-light scalar fields is well-known |l5j HH] and is very similar to the suppression due to 
free-streaming of light massive neutrinos. In this section, we review the origin of the effect and 
estimate the relevant scales and the amount of suppression for string theory axions. 

The axion field starts oscillating when the expansion rate drops below the axion mass at a 
time, 

H(rj m ) = m . (25) 

Before r} m the axion is frozen at a constant value, while afterwards it starts oscillating and its 
energy density redshifts away like ordinary CDM, 

( a m \ 3/2 r 

4>(v) = 4>o ~r~\ cosm / a(r])dr) . (26) 

\ a wJ J Vm 

A perturbation in the axion field, 5(p, with comoving momentum k satisfies 

5'<p + 2-5<p + k 2 5(j) + m 2 a 2 5(j) - 4# + 2m 2 a 2 0tf = , (27) 
a 

where \& is the perturbation of the gravitational potential (here we are working in Newtonian 
gauge and neglecting the anisotropic stress tensor, so that ^ uniquely characterizes scalar metric 
perturbations). At times later than r] m (and also such that the physical momentum of the mode 
is smaller than the mass, c.f. (33)) it is convenient to separate the oscillatory part of the axion 
field solution by writing 

5(f) = i ( e im V + e" im V*) , (28) 

where if) is a slowly varying function of space and time. Then the density perturbation in the 
axion field averaged over times longer than the period of oscillation is given by 

Pa \a m J 00 

and the axion field equation translates in the following equation for S a in the large mass limit 

k + % + (^-] S a = -k^ + 3^> + 3-vi> . (30) 
a \2ma j a 

This equation coincides with that for ordinary CDM apart from the last term on the lhs, which 
indicates the presence of a momentum dependent sound velocity, 

k 2 

c 2 s = . (31) 

4m 2 a 2 

This "quantum pressure" is a manifestation of the uncertainty principle for the axion particles and 
can be neglected if the physical momentum of a mode is smaller than the Jeans momentum \J Hm. 
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Assuming that the initial axion field is homogeneous, as implied by our choice of Hi n f ~ 0.1 GeV, 
S a satisfies adiabatic initial conditions and behaves just as CDM for scales larger than the Jeans 
scale. 

The evolution of a density perturbation of momentum k is ultimately determined by the time 
T) c of horizon crossing, 

~~ 7~~T = H( Vc ) , (32) 

and the time where the physical momentum becomes smaller than the axion mass and the mode 
is no longer relativistic, 

k , , 

= m . (33) 



a(r), 

Modes with comoving momentum 



, , m 1 ^ 4 ,/ m \ 1/2 
k < k m = ^ « 0.01 Mpc- 1 in _ 9S - (34) 



z 1 / 4 ' * V 4 x lO" 28 eV 



become non- relativistic while still being superhorizon ("long modes"), later the axion starts oscil- 
lating and finally the mode enters inside the horizon, 

Tj r < rjm < Vc ■ (35) 



Note, that in (34) we assumed that the axion is heavy enough, such that the beginning of oscilla- 



tions happens in the radiation dominated era r\ m -C rj eq . For these modes the additional term in 



equation (30) can be neglected and they behave like ordinary CDM. 

While the presence of an ultra-light axion component does not affect the CDM power spectrum 
is no longer true for modes with k > k m ("short modes"). These modes enter the 
horizon while still being relativistic, and they keep being relativistic at the moment 7] m , when the 
homogeneous axion field starts oscillating, so that the characteristic time-scales are ordered as 

Vc < Vm < f]r ■ (36) 



After r] r density perturbations in the axion component can be described by (30). The important 
difference with the long-mode case is that one cannot neglect the sound speed term in (30) until 
the even later time rjj, when the physical momentum of the mode becomes longer than the Jeans 
momentum (see Figure [2])), 

= v^M^ • (37) 

As a result any perturbation in such a mode decays until the moment r]j. The Jeans momentum 
redshifts with the same rate as the physical momentum of the mode during the radiation epoch, 



so that condition (37) gets satisfied for short modes only after matter-radiation equality. Since 
the momentum dependent sound velocity prevents density perturbations from growing, at any 
moment of time r] an axion behaves as a smooth homogeneous component of dark matter at 
physical momenta larger than ^jH{r])m. When an ultra-light axion constitutes only a subleading 
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Figure 2: Evolution of different physical momentum scales as a function of the scale factor. The 
horizontal long-dashed line corresponds to the axion mass. The short-dashed line shows Hubble. 
The dash-dotted line shows the evolution of the Jeans scale. Finally, three solid lines correspond 
to physical momenta of long and short modes and of the mode with k = k m . 



fraction of the dark matter the dominant consequence of such a behavior is a change in the growth 
rate of the dominant heavy CDM component after matter-radiation equality. Specifically, at this 
time the Poisson equation for the gravitational potential at scales shorter than the Jeans scale 
takes the following form 

-k^ = ^Ga 2 5p h « § ^^c?HH h , (38) 

where the subscript h refers to the dominant heavy component of the dark matter, and in the 
last equality we neglected the radiation and vacuum energy contributions to the energy density. 
Therefore the growth of the density perturbation Sh is described by the equation 

6 h + -5 h = 3 -^^a*H*5 h . (39) 

During matter dominance the growing solution of this equation evolves as 

5 h oc a p , (40) 

where 

-1 + J25 -24p°- 



4 5i l m 
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As a result, the presence of the ultra-light axion results in a step-like feature suppressing the dark 
matter power spectrum at small scales. The suppression starts at k m , and at shorter scales the 
suppression factor grows as 

S (zj) ~ ( i±^V 1 « i _ ^ i og I±*5 (42 ) 



l + zjj 5fl m l + zj 

where z eq « 3200 is the redshift at matter-radiation equality, and zj is the redshift corresponding 
to the moment rjj when the physical momentum of the mode becomes longer than the Jeans scale 
kj = y/Hm, 

1 + *J ~ fc4 ° • (43) 

For modes with momenta higher than kj at the redshift of observation, z , the effect saturates; the 
power spectrum for such modes is suppressed relative to the usual ACDM case by a factor S(z ). 
Consquently, the presence of a subdominant ultra-light axion component exhibits itself as a step- 
like feature in the power spectrum as shown in Figure |3j The width of the step is almost an order 
of magnitude in the comoving momenta — the suppression shows up at k m ~ (mHo) 1 ' 2 (Q m /z eq ) 1 ' i 
and saturates at around kj ~ (mi/o) 1 ^ 2 (^m) 1//4 - As one might expect the magnitude of the effect 
is controlled by the fraction Q a /Q m of the axion density relative to the total CDM density, however 
there is also an additional logarithmic enhancement, 

log^^«8-log(l + , (44) 

due to the accumulation of the effect over a long period of time. This discussion assumes z Q <C z eq , 
at z a ~ z eq the suppression factor is simply ~ Q a /Q m . 

It is straightforward to estimate the fractional abundance Q a /Q m : An axion starts to oscillate 
when H ~ m, and at this time the total energy density in the Universe is equal to 

p tot ~ 3M 2 Pl ^ , (45) 

J a 

while the axion density at the beginning of oscillations is of order A£. As a result one finds that 
the axion fraction is of order 



HOi) ■ • , (46) 

n m 1 l> 3M|, 1 + z eq ' K J 



where 



1/2 1/4 / \ 1/2 

m i z eq ( m \ 



H 1/2 Q 1/4 9 \4xl0- 28 eV 



(47) 



is the redshift at H = m and P{0i) is a statistical factor depending on the initial axion angle 



Qi — 0/ 'fa- We see that the size of the steps, which is the product of the fraction (46) and 



a logarithmic factor (44), depends on three parameters f a and z m so we now discuss the 
dependence and typical values we expect for each of them in turn. 

We illustrate the dependence on the initial axion angle in Figure [4] (c.f. [17]). On the left panel 



we present solutions for an axion field with different initial conditions in a radiation dominated 
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Figure 3: Suppression of the power spectrum observed today as a function of the co moving 
monentum. 5u has been normalized to the value at large scales and we have assumed 5 s - = 0.01. 



Universe as a function of a scale factor normalized to one at H = m. We see that increasing the 
initial value of the axion gives rise to two effects. First, the resulting axion amplitude is higher 
just because the initial condition is higher. Second, the axion potential becomes flatter close to 
9 = 7r, and as a result oscillations start later also resulting in the increase of the axion abundance. 
We illustrate this effect on the right hand panel of Figure [4j where we present a numerical result 
for the enhancement factor P(9) as a function of the probability (tt — 9) /it to have an initial axion 
misalignment higher than 9. 

The crucial microscopic parameter that determines a size of the effect is the ratio /f/(3M|J. 
As discussed in Section 1.2 , it is widely believed that this ratio is necessarily smaller than one, and 
the relation f a /Mpi < 1/S should hold, where for the ultralight axions that we consider here the 
instanton action should be quite big, S > 200. This favors rather small values /^/(SMj,;) ~ 10~ 5 . 
On the other hand, these estimates have an order one uncertainty in them, so that values as high 
as fa/(3Mp t ) ~ 10~ 4 may be not unreasonable. 

Finally, the third factor in (46) — the redshift ratio — is determined by the axion mass. Different 
types of observations will be sensitive to different values of masses. For instance, CMB measure- 
ments are sensitive for the CDM power at scales crossing the horizon around matter-radiation 
equality. From the above discussion we see that there is highly unlikely to be an effect at the level 
significantly higher than 10 -3 at these scales (note also, that in this case z Q ~ z eq so one looses 
the enhancement of (44)), so the CMB is not a good probe of this effect. 

At the shorter scales the limits on the CDM power spectrum can be obtained from galaxy 
surveys and Ly-a forest data. The step-like feature in the power spectrum due to ultra-light 
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Figure 4: A time evolution of the axion field for different initial conditions (left panel) and an 
enhancement factor P(6) for the axion abundance as a function of a probability for different initial 
axion values (right panel). 



axion is very similar to the suppression of the small scale power due to a free-streeming neutrino 
component (c.f. [46J), so to estimate the sensitivity of the currently available data we can translate 
the bounds on the warm-plus-cold dark matter models involving sterile neutrinos (see jlH] for the 
up-to-date analysis). The bounds are dominated by the SDSS Ly-a data, which are sensitive 
for the comoving momenta in the range k ~ 0.1 -v- 10 Mpc -1 and correspond to the observation 
redshift z ~ 2 4. For axion masses close to the lower end of the above momentum range this 
analysis allows admixture of an axion component at the level Q a /Q m ~ 0.1, while for masses 
corresponding to k > 4 Mpc -1 an order one fraction of an ultralight axion component is allowed. 
Taking into account that matter-radiation equality corresponds to k eq ~ O.OlMpc -1 , and that the 
characteristic comoving momentum k m scales linearly with z m we conclude from (46) that these 
bounds are still not good enough to probe the most interesting region /^/(3Mp ; ) ~ 10~ 5 . 

On the other hand, the Baryon Oscillation Spectroscopic Survey (BOSS — which is a part of 
SDSS III [35]) will have a sensitivity at the level of few percent to the CDM power spectrum at 
comoving momenta k ~ O.lMpc -1 , which correspond to axion masses around 4 x 10~ 26 eV (cf. 
(34)). Given that matter-radiation equality corresponds to k eq ~ O.OlMpc -1 the redshift ratio 
gives a factor ~ 20 at these scales (as follows from Figure lik) the axion starts oscillating at a 



redshift at least half that of z m ). For these observations z Q ~ 1 so the log enhancement (44) is 
maximal. Therefore, with a reasonable statistical factor, P{9) ~ 20, BOSS will be able to observe 
the effect of axions down to fl/C&Mpt) ~ 10~ 5 . 

In the more distant future, a high precision measurement of axion steps from 21 cm line 
measurements is possible. Indeed, observations of the 21 cm line have a chance to probe the power 
spectrum in the range k ~ 10~ 2 4- 10 3 Mpc^ 1 These observations correspond to z Q ~ 30 -j- 200, 
so one looses slightly the enhancement factor (44). However, this will be easily compensated by 
the increase of z m and by the high accuracy that can be achieved by 21 cm measurements. Also, 
unlike for the Ly-a forest, this has an advantage of probing the power spectrum in the linear 
regime. To appreciate the precision one may hope to achieve, it is enough to note that the 21 cm 
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line will provide a total of N21 ~ 10 16 independent modes as opposed to N cm b ~ 10 7 modes for 
the CMB. 

Note, that axions with lower masses tend to produce more pronounced steps, so that at the 
early stages when only the low k data is available, rather than observing an individual step one 
may discover an overall running of the spectral index due to a joint effect of several axions with 
different masses. 



As follows from (46) an important scale is k ~ 10 4 : 5 k eq — where axions with the corresponding 



masses are expected to constitute an order one fraction of the CDM for a typical initial misalign- 
ment. Axions with heavier masses would overdose the Universe for a typical misalignment, so 
their initial amplitudes are forced to be fine-tuned to small values in our patch of the Multiverse. 
On the other hand, there is no reason for initial values of the lighter axions to be fine-tuned, so 
we expect that axions with masses close to the anthropic boundary, i. e. of order 

P{9)' V fit ) HO) 2 \ 10* I 1 ' 

to constitute an order one fraction of the CDM. The 21 cm line measurements can probe part 
of the anthropic regime, as they are sensitive to masses up to 3 x 10~ 18 eV. The level of the 
expected signal is uncertain though, and may depend on the probability measure that determines 
the fine-tuning of the initial axion amplitudes. 



The value of the axion mass in (48) is also interesting as it has been previously suggested 
that the problem with excessive small-scale structures predicted in the vanilla ACDM model may 
be solved if dark matter is composed of a Bose-Einstein condensate of ultra-light particles with 
masses of ~ 10~ 22 eV. Of course, it may be that the problem will be resolved within conventional 
ACDM after more precise iV-body simulations including baryons become available. Still, it is 
intriguing that a similar mass scale appeared in the present context as the boundary between 
anthropic and non-anthropic axions, so that one expects an order one fraction of CDM to be 
composed of axions in this mass range (recall that the statistical factor of P(0) ~ 10 has a 
probability 0.2 for a single axion). Note, that this connection, if true, favors high /2/(3Mpj) 
ratios around 10~ 4 . On a positive note this increases the likelihood of seeing light axions with 
BOSS, but it forces a reconsideration of whether such a high ratio is achievable for ultra-light 
axions in explicit string theory constructions. 



2.4 Extraction of Black Hole Rotational Energy by Axions 

Ultralight axions with Compton wavelengths comparable to the size of astrophysical black holes 
can be sought for in observations of rapidly spinning black holes. The key reason for this is the 
Penrose process, which opens the channel for the black hole spindown. Indeed, the spinning black 
holes have the so-called ergoregion, a region outside of the event horizon and therefore accessible 
to external observers, inside which the inertial frame dragging due to the black hole rotation is 
so fast that probes built of normal causal matter, which can never move faster than light, cannot 
remain at rest relative to an observer far away. This opens the possibility for energy and spin 
loss by the black hole |1HJ \5U\, which happens when a particle falls into the ergoregion such that 
it co-rotates with the black hole. If after diving into the ergoregion the projectile splinters up 
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into two fragments, and one falls into the black hole while the other recoils away and out of the 
ergoregion, the escaping fragment can take out more energy and spin than the original projectile. 

To see this in a bit more detail, consider the geometry of the spinning black hole given by the 
Kerr solution in Boyer-Lindquist coordinates (51] 



, 2 .„ 2R q r. l2 4iLar sin 2 9 , ,, E l2 „ in2 (r 2 + a 2 ) 2 - a 2 A sin 2 9 . , 

ds 2 = -(1 ^~)dt 2 9 — dtdcf) + —dr 2 + Zd9 2 + J —— sin 2 

Zj Zj Za Zj 



S = r 2 + a 2 cos 2 #, A = r 2 - 2R g r + a 2 , a = j^, R g = G N M , (49) 



and where M and J are black hole's mass and spin, respectively. The roots of £ = 2R g r define the 
ergosphere, whose interior is the ergoregion, while the roots of A = are the black hole horizons, 
and the outer one, relevant for our discussion, resides at r + = R g +y/R 2 — a 2 . The spatial location 
of the outer horizon never extends past the ergosphere, reaching it only at the poles. The time at 
asymptotic infinity is measured by t, and represents the clock reading of an inertial Minkowskian 
observer far from the black hole. The time evolution of a physical system is therefore described by 
the flow generated by the vector field 7i = d t . Because this vector field is also a Killing vector of 



the geometry (49), if we consider an inertial particle of mass fi in this background, the product of 
its 4-momentum p M = dx^/dr (where dr 2 = —ds 2 ) and of the Hamiltonian is conserved, yielding 
the energy integral of motion, 

E = -Hp . (50) 



Similarly, since (49) is manifestly axially symmetric, and the rotations around the z-axis are 



generated by the vector field J = d^, the product 

L = Jp (51) 

will be the conserved particle's angular momentum in the z-direction. 

The crucial property of the ergoregion is that particles moving there may have a negative 



energy. Indeed, inside the ergoregion the (it)-component of the Kerr metric (49) is negative, so 



that the Killing vector 7i is space-like. The 4- velocity of any physical observer should be time-like, 



implying that all observers experience rotation inside the ergosphere, d^jdr > 0. The energy (50), 
being the product of a time-like and a space-like vectors, is not sign-definite, and takes negative 
values for some of the observers. 

This observation is the key to the Penrose process of black hole energy loss. Imagine a projectile 
falling toward the black hole, initially coming in from afar where 7i is time-like. Then its total 
conserved energy must be positive. Allow it to dive beneath the ergosphere, but aim it such that 
it misses the horizon. Then design a timer on the projectile to set off a fragmentation process 
that will break the projectile into two parts, but so that one fragment flies off into the black 
hole along a trajectory on which the conserved energy is negative, £ < 0. Since the total 4- 
momentum conservation governs the fragment dynamics, the energy of the escaping fragment will 
be E = E — £ = E + \£\, exceeding the energy of the initial projectile. Its gain comes at black 
hole's expense: the conservation of the total energy of the system implies that 

8M = -{E-E) = £ < 0. (52) 
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As a result of this process the rotation of the black holes spins down. Indeed, the vector field 
Q = TC + uj + J is a Killing vector that becomes null on the horizon. Here ui + = a/(2R g r + ) is the 
angular velocity of the horizon. For later purposes, we note that as a function of the spin a and 
the gravitational radius R g , it is given by 

= 1 a / R 9 (53) 

+ 2R g 1 + VI - (a/R g y 

Note, that in the extremal limit a/R g — 1 it saturates at R g u> + (ma,x) = 0.5. Now, the product of 
the 4-momentum V of the infalling fragment (which is a future directed time-like vector) and of 



the future-directed null vector Q must obey Q ■ V < 0. Using the definition Q and Eqs. (50) and 



(51), this yields for the spin of the infalling particle C < £/lu + . Consequently, by total angular 



momentum conservation, there will be a change of the black hole's spin as well, 

5 J = -(L-L) = £ < — < 0. (54) 



Eqs. (52) and (54) show that for Penrose fragments, both the energy £ and the spin C are negative. 



Combining these two equations we find that the black hole hair change by 

SJ<^, (55) 

UJ + 

which in fact is an equation that hides in itself the second law of black hole thermodynamics. 



Indeed, one can compute the horizon area of the spinning black hole (49) to find 



such that 



'l-a 2 /R 2 p 

A = lGnB* , (56) 



^/Rj^\u + 



By Eq. (55) this yields 



SAX), (57) 

precisely the area law. In principle, it is possible to fine-tune the fragment trajectories in the 
Penrose process such that S A = 0, and so mine out black hole energy stored in the spin in the 
amount of 

M ( 



5£ = ^\^-^l + y/l-a*/B* \ , (58) 

which means the more energy will come out if there was more spin to begin with. The extraction 
is maximized for (near) extremal black holes, where it reaches 29% of the total mass. Note that 
the energy stored in spin is a very sensitive function of a/R g : for example, for a/R g ~ 0.8, the 
energy is already significantly smaller, 5£/M ~ 0.106, whereas for a slowly spinning black hole 
with a/Rg ~ 0.3, it is 5£/M ~ 0.0116. This energy is essentially analogous to the difference of 
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the total 'mc 2 ' energy and rest energy rriQC 2 in special relativity, and shows that the black holes 
with spins a/R g > 0.9 are very relativistic, and so fast-spinning, those with a/R g ~ 0.8 or so are 
analogous to 'warm' particles with kinetic energies comparable to the rest mass, while those with 
smaller spins, e.g. a/R g ~ 0.3 are already 'slow'. 

The above discussion was presented in terms of point particles, while we are interested in 
applying the Penrose process to scalar fields — axions. One way to see that one can use them to 



extract energy from the black hole is to reverse the logic and note that the inequality (55) is a 



direct consequence of the area law (57), which is nothing but the second law of thermodynamics. 
Consequently the inequality (55) should hold independently of what carries the energy and spin 
into the black hole. Applied to the incoming scalar wave of the form \l/ = e - l ^ t + tm( t> j (j- ^ with a 
frequency in the range 

< u < mw + (59) 

this inequality implies that both energy and spin transfer from such a wave into a black hole 
should be negative [S2J G21 IM]- A direct way to see this [SI] is to consider a conserved energy 
current for the scalar field \l/ given by 

p» = -t^u v = -dfld& ■ 

Let us now consider a space-time region between two slices of constant time t. Conservation of the 
current P M implies that the time-averaged energy flux at the infinity is equal to the time-averaged 
energy flux through the black hole horizon. The latter is equal to 

W> = + u+d+tydt*) = lo{u - muj + )f 2 (60) 



and is indeed negative in the frequency range satisfying the superradiance condition (59), so that 
by scattering such waves off the black hole one extracts the energy. Based on this observation 
Press and Teukolsky [23 EE] designed a "black hole bomb" — a spinning black hole surrounded 
by a spherical mirror. Such a device exhibits an exponential classical instability — being confined 
by the mirror a small initial scalar field perturbation inside the shell with a frequency in the 



superradiant range (59) experiences a repeated series of amplifications by scattering off the black 
hole, until it extracts all of the black hole's spin (or the mirror blows up). 

An extremely interesting observation made already in [57] is that the Nature itself may provide 
such a mirror if the field has a non-zero mass. Indeed, unlike massless particles, massive ones can 
rotate on stable orbits around the black hole just like planets around the Sun. Consequently, for 
a massive scalar field there should be a set of bound states in the Kerr background corresponding 
to wave-packets rotating along these stable orbits. However, unlike for the point particles, such 
a wave packet will always have a tail going into the ergosphere region as well. If the wave- 



packet contains frequencies in the superradiant range (59) they will be continuously amplified 
and the amplitude of the field will be growing exponentially. In other words, one expects to find 
an exponentially growing bound states in the spectrum of scalar field perturbations in the Kerr 
background. This intuition was proven to be correct [58l [59] . Let us see how unstable modes arise 
at a more technical level. Rewriting (□ — fj 2 )^ = in the Kerr background (49) and using the 
separation of variables *Sf = R(r)Q(8)e irn ^e~ llJJt , one finds that is an oblate spheroidal harmonic, 
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and that R must satisfy 



A-^- (A^j?) + (a 2 m 2 - AR g ramu + (r 2 + a 2 ) 2 uj 2 - /iVa) R = (X ml + cA 2 ) AR , (61) 

where X m i are oblate spheroidal eigenvalues, which depend on m and I but in general cannot 
be written analytically in terms of them. This equation can be cast in a more useful form as a 
Schrodinger equation, by defining the tortoise coordinate r* according to dr* = (r 2 + a 2 )dr/A, 
and rescaling the wave function to u = \Jr 2 + a 2 R. The resulting equation is 



+ (^ 2 -VH) M = O, (62) 



with the potential 



A/i 2 ^ AR g ramu — a 2 m 2 + A (A m « + (cj 2 — [i 2 )a 2 ) 
r 2 + a 2 (r 2 + a 2 ) 2 

A(3r 2 - AR g r + a 2 ) 3A 2 r 2 
+ (r 2 + a 2 ) 3 (r 2 + a 2 ) 4 ' ^ 

The potential approaches — > Voo = fi 2 far from the hole, when r* — * oo, rising out of a potential 
well where V < jj 2 (see Figure [H]). Note, however, that it rises towards its asymptotic value at 
infinity already as soon as r > and so the effective spatial extent of the potential well is 
really 8r ~ 1///. This "mass barrier" plays a role of the mirror which reflects the runaway Penrose 
fragment back toward the black hole, enabling it to undergo repeated Penrose scatterings and 
gain more energy and spin from the black hole. On the other side of the well is the centrifugal 
barrier, peaking at r* ~ R g . Beyond the centrifugal barrier is the ergoregion, where the potential 
asymptotes to V — > V + — 2muj + uj — m 2 u; 2 as r — >• r + (equivalently, as r* — > — oo). 



The crucial property of the equation (62) is that it is not a selfadjoint eigenvalue problem- 



the potential (63) depends on u in an essential way when a ^ (equivalently co + 0). This 



makes possible the existence of bound states with complex frequencies. In particular, near the 



horizon the terms in the parenthesis in (62) combine into uo 2 — V + = (a? — mw + ) 2 . As a result, the 



asymptotic form of the wave function near the horizon is exp(±ifc + r*), where 

k. = uj — ma), . 



The requirement of regularity at the horizon singles out one of these waves, exp(— ik + r*) 
The appropriate boundary condition at infinity, given that we are looking for bound states, is the 
exponential decay of u. Altogether, the boundary conditions are 

^ _^ < S>ml(0) e - ik+ r*+im4>-iut ^ for r * _^ _ QO ^ (g 4 ) 



^2R : 



'3^ 



^ _^ Qm;(^) c -^/^-u)2 r *+im<f>-iu)t 



for r* — > +oo . (65) 
r 
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Figure 5: The effective potential of Eq. (63). Depicted are the ergoregion, to the left, with the 
horizon at r* — > — oo, the centrifugal barrier (whose height depends on the angular momentum 
of a mode), the potential well to the right of it, and the asymptotic mass barrier which plays the 
role of the mirror that reflects the escaping Penrose fragment back. The relevant modes will be 
the states bound in the potential, and leaking through the barrier towards the horizon. 



We see, that bound states with Re k + > satisfy ingoing boundary condition at the horizon, 
i.e. being unable to escape at the infinity, they still can be sucked into the black hole. On the 



other hand, for Re k + < 0, i.e. for frequencies in the superradiant interval (59), the boundary 



condition (64) describes an outgoing flux of particles from the black hole. Of course, this is 
in a perfect agreement with the previous derivation demonstrating that superradiant modes get 
amplified as a result of scattering off the black hole. This implies that if the real part of a 
frequency is smaller than mw + , the imaginary part should be positive, indicating the presence of 
an exponential classical instability. In the Appendix we worked out the details of how this happens 
in a simple toy potential. The quantitative expressions for the imaginary parts in the actual Kerr 
background were caculated explicitly in the limit R g jj, ^> 1 in [58] and in the limit R g fi -C 1 [59] . 
and the numerical calculation in a general case was performed recently in [60]. With the notation 
uj = uj t + i/r, the explicit calculations find 

T = io 7 e L84i? ^ R g ) f or R gfJi > i an d a = \ ; (66) 

•9 



t = 24 (J- ) ( /?„,/, fur /?,,,, «1. (67) 



For a < 1 the instability scale at large masses has the same qualitative form (66 ), but the coefficient 
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1.84 in the exponent grows at smaller a. The following comments are in order regarding the fastest 
instability channel. First, the instability 



benefits from as large a as possible, being the fastest for near extremal black holes with 
a/R g = 1. 



Since for the unstable modes tu r < mu + , and uj + < l/(2R g ) (see Eq. (53)), we find that u r < 
m/{2R g ). Now, the rate of the instability is controlled by the tunneling of such modes through 
the centrifugal barrier, which is higher at larger I. Hence, among the available modes the fastest 
growing one will be the one which tunnels most easily, i.e. 

• has the smallest total angular momentum I; 

• has the largest projection on the black hole's axis of rotation, m — I. 

The first condition implies that it has the lowest possible barrier, and the second ensures that it 
is the highest level in the well, probing the thinnest available section of the barrier. Now, the 
specific frequencies that can be used depend on the value of /i as well. If R g fi ^> 1, then fi ^> uj + , 
and the fastest instability is due to the mode with m ~ I ~ resulting in the exponential 

suppression due to the fact that the potential well is very narrow and the centrifugal barrier high 
and thick. In the opposite limit, R g fi <C 1, the fastest instability is in the sector m — I — 1, which 
is now available, and is suppressed mostly due to the fact that the fastest mode is spread out 
through the potential well becoming wider near its top. This results in the power law suppression 
of the instability rate. 

Numerical results of [60J show that the instability peaks around R g n ~ 0.42 (for a/R g ~ 0.99), 
where 

r sr ~ 0.6 x 10% , (68) 

In this regime the instability is due to the mildly non-relativistic (the real frequency is u r ~ 0.98/u) 
I = m = 1 level. 



2.5 Axionic Sirens and Precision Black Hole Physics 

Black holes are believed to be abundant in our Universe; in particular, there is a 4 x 10 6 solar mass 
(Mq) black hole at the center of our Galaxy. Current black hole candidates are primarily found 
in X-ray binaries and Active Galactic Nuclei (AGN) and clustered in two mass ranges, 3M - 
3OM for stellar mass black holes, and 1O 6 M — 1O 9 M for supermassive black holes [HI]. There is 
currently an evidence also for intermediate mass black holes [62] as well as for black holes heavier 
than 1O 1O M [53]. 

In the past few years, a lot of progress has been made in measuring the properties of these 
black holes including their spin. There is evidence for rotating stellar mass black holes as well 
as a rotating black hole candidate at 1O 7 M \§T\. In the next decade X-ray observations, in 
combination with a better understanding of back hole environments, will solidify and extend 
black hole spin measurements. In the not-so-near future gravity wave observatories such as LISA 
and AGIS [62] will provide even more precise probes of black hole properties. The natural question 
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arises what implications the axion superradiance may have for the fate of astrophysical black holes 
and whether observations of astrophysical black holes may provide evidence for ultra-light axions. 
To address this question the analysis of section |2.4 has to be extended in several directions. 



First, one may worry that the instability may be significantly inhibited or even totally disappear 
in a realistic astrophysical enviroment. Apart from the exponentially unstable superradiant modes 
Kerr metric also supports bound states in the non-superradiant regime, whose frequencies have 
negative imaginary parts. Physically, these bound states get damped rather than amplified by the 
black hole and eventually get sucked behind the horizon. Some of these modes may have imaginary 
parts much larger than the maximum value r^, 1 ~ lO -7 /?" 1 for the rate of the superradiance 
instability [60J - Consequently, one may worry that perturbations to the Kerr metric which are 
always present in a realistic astrophysical situation may lead to the mixing between superradiant 
and non-superradiant modes and damp the instability. 

Second, if the superradiant instability survives these perturbations and developes, at some 
point the linearized approximation breaks down, and the backreaction of the axion cloud has to 
be taken into account to deduce the observational consequences. 

2.5.1 Superradiance in Realistic Enviroments 

Let us start with checking that the superradiant instability persists also in a realistic astrophysical 
enviroment. We will see that the black hole vicinity is a very clean astrophysical enviroment — as 
exemplified by no-hair theorems — the black hole itself cleans up the space around its horizon. 

One source of perturbations on the Kerr metric comes from the presence of accreting matter. 
A reasonable estimate (in fact, an upper bound in most cases) for the accretion rate is provided 
by the Eddington limit [65] — the accretion rate, such that the radiation pressure on free electrons 
balances gravity. In the Eddington regime the black hole mass M grows according to 

M 

M = —— , (69) 

Taccr 

where the Eddington accretion time T accr is equal to 

racer = , !T mpson ~ 4 x 10 8 years. (70) 
47rGArm proton 

From this we can estimate the amount of matter in the vicinity of the black hole horizon as 

SM ~ M^- « 4 x 1(T 22 M ( —) . (71) 

The shifts of the imaginary parts one can expect from such a perturbation are of order R~ X 8M jM , 
which is comfortably smaller than the positive imaginary part due to superradiance, if the axion 
mass is not too far from the optimal regime. 

Another source of potentially dangerous distortions is a tidal force due to a companion star 
rotating around the black hole. For instance, stellar mass black holes are observed in X-ray 
binaries, so such a star is always present. The physical (gauge invariant) part of the metric 
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perturbation due to tidal forces caused by a companion star with a gravitational radius r g at a 
distance L is of order 

where the Riemann curvature created by the companion star is of order Rn\ up ~ r g/L 3 and 
Aa; ~ R g is the size of the axion cloud. Consequently, the tidal force correction to the imaginary 
part of the superradiant modes is at the level 

fa ~ ^tf-i . (72) 

To get an idea of how big this correction is let us estimate it for the X-ray binary LMC X-l, 
supposedly harboring a 1OM black hole with spin a/R g ~ 0.91 [66] (the motivation for this 



choice is that, as we discuss in Section 2.5.3, at the moment LMC X-l provides the best reach for 
the QCD axion). The companion star in the LMC X-l is quite heavy — around 30 M |67J, and 
its orbital period is ~ 3.9 days, which corresponds to the distance L ~ 2.5 x 10 12 cm. Altogether, 
this gives rise to a tiny frequency shift of order ~ lO -18 /?" 1 , which is not dangerous for the 
superradiance. 

As an example of a galactic black hole let us consider Sgr A* — a ~ 4 x 1O 6 M black hole at 
the center of the Milky Way. For a number of years Keplerian orbits of around 30 stars around 
Sgr A* were monitored [681 . with the fastest period being equal to 15 yr, which corresponds to the 



distance of order 10 4 R g , again being safe for superradiance instability. It follows from (72) that 
for relatively light galaxies, M < 1O 7 M , a compact stellar size object (such as a neutron star, or 
a black hole) within ~ 10R g from the horizon may be dangerous for the superradiance instability. 
However, such occurances are likely to be quite rare, as follows from the estimates for the rate of 
extremal mass ratio inspirals (EMRI), ~ few x 100 Gyr~ x [69] and are likely not to last long, as 



the objects get swallowed by the black hole. Given that the time (68) to build up the axion cloud 
is only ~ 10 4 yr even for the largest (1O 1O M ) black holes, EMRI's effect on the superradiance is 
also negligible. 

Another potential source for the superradiance shutdown, could be the presence of a magnetic 
field around the black hole if the axion has a coupling to EM 

Cn ' ^ Xp F, u F Xp , (73) 



47T/, 



where C is a constant that is ^ in SU(5) GUTs. As the axion instability develops around the 
black hole, there will be conversion of axions to photons due to the presence of a strong magnetic 
field and energy will be carried away from the axion field. A strong enough magnetic field may 
lead to a total dissipation of the axion cloud. 



In the presence of non-zero magnetic B and axion fields the coupling (73) gives rise to the 
source term in the Maxwell's equations 

= ^re^F, p d x 4> ~ ^B<f> , 

2vr/ a W Tlf a 
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where we took into account that the axion field oscillates with a frequency of order m. This source 
produces photons with an energy density of order 



that take away the energy from the axion field with a characteristic time scale 



T ~ 

'7 



f \ 2 
71" fa \ /X 



Cat J W 2 

In order for the superradiance instability to be effective r 7 should be longer than the characteristic 



superradiance time r sr (68) or, equivalently 



B<2- 1(T 4 ^A 2 , (74) 

where A = y/mf a , assuming ~ (CrjvM) -1 . 

For the QCD axion the bound (74) corresponds to a magnetic field of order 5 • 10 16 G, which 



is much larger as compared to what the black hole accretion disk can support (Eq. (|75|), [65 



M 



-1/2 



B ~ 4xl0G UJ ■ (75) 

Hence the superradiance instability is definitely present, and the magnetic field does not affect the 



range of / a 's probed. Similarly, (75) implies that (74) is satisfied for supermassive galactic black 
holes as well for f a ~ M GUT . 

However, magnetic fields may shut down the superradiance for f a <C Mqut- An example is 
the ultralight axion which can affect supernovae luminosities by photon-axion mixing [70J, and 
which operates at the mass scale \i ~ few x 10~ 16 eV with the effective coupling Ca/(47if a ) ~ 
(4 x 10 11 GeV)~ 1 . This field could be in the regime of the fastest instability of the smallest 
supermassive black holes on the record, with M ~ 1O 5 M , which are able to support magnetic 



fields B ~ 10 6 G, while the critical S-field value for them, by Eq. (74), is ~ 10 2 G. Thus the 
instability may be turned off in this case. 

2.5.2 The Fate of the Axionic Instability 

Let us discuss now what happens with superradiance at late stages of the instability development. 
As always with a linear instability, at a certain point the backreaction of the axion cloud has 
to be taken into account and the problem becomes non-linear. Naively, one may think that at 
this stage the problem becomes very complicated and hardly tractable. However, the important 
property of the superradiant instability is that it is always very slow — the instability time scale 
is at least seven orders of magnitude longer than a natural dynamical time-scale of the system 



R g , see (68). As a result, second order effects have a chance to compete with the superradiance 
instability when the axion cloud is still a small perturbation to the black hole, and the whole 
process can be under control. Still, as we will see, a black hole surrounded by an axionic cloud is 
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Figure 6: The Carnot Cycle of the Axionic Instability: The black hole "feeds" the su- 
perradiant state forming an axion Bose-Einstein condensate. Axions from that state quantum- 
mechanically transition through graviton emission to a lower-energy non-superradiant state. The 
non-superradiant state decays into the black hole. Accreting matter around the black hole replen- 
ishes the rotational energy lost to gravitons and sustains this cycle. 



an extremely rich dynamical system, and our purpose here is just to provide a basic intuition of 
what one can expect from it. A dedicated quantitative analysis of different regimes will appear in 
a separate publication [71] . 

A rotating black hole surrounded by an axionic cloud is essentially a huge quantum mechan- 
ical system similar to an atom, with a crucial difference that particles populating its levels are 
bosons, rather than fermions. As a result, some of the levels become highly populated. More- 
over, a nucleus — the central black hole — continuously creates particles on some of the levels and 
destroys them on the others with the rate proportional to the occupation number of the level. A 
model independent source of back-reaction is related to the possibility for an axion in the growing 
superradiant cloud to emit a graviton and to jump onto a non-superradiant level, from where it 
eventually gets sucked inside the black hole horizon. The whole system is fueled by the inflow 
of accreting matter, and, if the accretion is efficient enough, it acts as a giant gravitational wave 
siren, see Figure |6j 

To gain some intuition about different regimes of the axionic siren it is instructive to write a 
simplified set of kinetic equations, describing the evolution of the black hole, accretion disc and 
axionic cloud. Namely, let us characterize the siren by the total number of axions on the superra- 
diant levels N + , non-superradiant levels AL, and the black hole spin J = jh- The evolution of the 
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occupation numbers N+, AL can be approximated by the following couple of kinetic equations, 



d ^f = r^ 1 N + -T G l v N + (N_ + l) (76) 
^- = -rZ 1 ^ + T£yN+(N- + 1) (77) 

Here, 

r± = t±R 9 (78) 
are the superradiance and dumping time, and the graviton emission time tqw is 

r GW ~ M 2 Pl R 3 g , (79) 

assuming the axion mass is close to the optimal value, fi ~ R~ 1 . This two-level approximation 
becomes literally true in the low mass limit fi < -R" 1 , when the axion cloud is dominated by 2p 



states. The pair of kinetic equations (76), (77) has a stationary point, at which 



N ± = T^K ~ e^M 2 Pl R 2 g = e T MR„ . (80) 

As we expected, in this stationary regime the total mass of the axion cloud M a ~ n{N + + AL) is 
small compared to the black hole mass. In order to establish this regime, the accretion of matter 
into the black hole should be efficient enough, otherwise the black hole would simply loose its spin 



and the siren would not start. To see when the stationary regime (80) can be established for a 
rapidly rotating black hole, let's consider the following kinetic equation describing the evolution 
of the black hole spin, 

dJ * /o"i \ 

~ir = T accrJ - Jgw , (81) 

where we assumed the Eddington regime for the spin accretion, so that the accretion time is given 



by (70). The spindown rate due to gravitational wave emission J GW at the stationary point (80) 
is determined by the number of gravitons emitted per unit time (each graviton carries away two 
units of the angular momentum) 

Jgw « 2t^N+N- ~ MR g € -^ = e^r" 1 . (82) 

Kg 

Consequently, the axionic siren can operate in a stationary regime provided 

T sr > e_T accr . (83) 

Note, that this condition is different from the naive one r sr > r accr , the reason being that the 
axion cloud is small, so that less accretion is needed to support it. In the opposite regime, 
r sr < €-T accr the accretion is not efficient enough to support the siren, and the black hole spins 
down. Consequently, in the presence of light axions we expect to see gaps in the spectrum of 
rotating black holes at the masses close to the optimal. We illustrated this effect in Figures [7], 
|H In Figure ft] we show the regions in the black hole mass/axion mass parameter space, where 
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Figure 7: The parameter space in the plane of axion and black hole masses where the superradiance 
leads to the black hole spindown assuming Eddington accretion (dark shaded region), and no 
accretion, i.e. the superradiance time scale is required to be faster than the age of the Universe 
(light shaded area). In both cases the upper bound has been calculated from (66) and the lower 



bound from (67). 



the condition (83) cannot be satisfied for a maximally spinning black hole and, consequently, 



black holes cannot sustain maximal rotation. The dark shaded region arises for Eddington limited 



accretion, see (70), while the light shaded region for an accretion time of order the age of the 
Universe, i.e. currently non-accreting black holes. In this plot we assumed that e_ = e + and that 
r sr given by (66 ) and (67) for the heavy and light axion mass regime, respectively. This plot should 
be considered only as an indicative one and the more refined analysis will be presented in [TTJ . In 



particular, the estimate (79) for the graviton emission time should be considered as a lower bound, 



as it does not take into account non-relativstic suppression. Taking this suppression into account 
will result in the stronger bounds on the axion masses than the conservative estimates presented 
here. We checked, however, that setting e_ = 1 in (83) affects the allowed axion masses in the 



heavy region only by a factor of order two (while the lower boundary of the spindown region in 
Figure [7] is more sensitive to such a change). 

In Figure [8] we illustrated how the black hole Regge plane (the parameter space of black hole 
masses and spins) may look like after precision measurements of spins and masses for many black 
holes will become available. Namely, in the presence of two axions with masses 1.7 x 10~ n eV 
and 3 x 10~ 17 eV black holes will populate only the shaded region in this plot (there could be 
rare exceptions due to black holes experiencing a short period of superEddington accretion). This 
figure should be considered as indicative, and is meant to illustrate that the dip in the spectrum 
of rotating black holes becomes narrower and less pronounced with increasing black hole mass. 
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Figure 8: The maximum allowed spin for a black hole as a function of its mass assuming there 
are two axions with mass m a \ and m a 2 corresponding roughly to black hole masses of 2M and 
1O 6 M . This plot has been created using (66) (and the dependence of the superradiance rate on 
the black hole spin for heavy masses from |58j) and (67) and is indicative. 



Heavy enough black holes may satisfy the condition (83) and operate in the siren regime. The 



graviton flux (82) corresponds to the gravitational wave signal at the Earth of the strength 
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where L is the distance to the siren. For e± = 10 this translates into 



3 x 10 
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which is above LISA sensitivity. 



We should stress, that the stationary point (80) is the simplest possible regime for the axionic 



siren. However, the siren is a very rich and complicated dynamical system that may exhibit other 
even more colorful periods during its life-time. For instance, if we deviate from the stationary 



point (80) of equations (76) and (77), we find cyclic solutions with the axion population oscillating 
between the superradiant and the non-superradiant level. These solutions would give rise to a 
gravitational pulsar with a period of order the superradiance time (10 years for a 1O 7 M black 
hole), with a maximum amount of radiation during the periods when the population of both levels 
are comparable. A cyclicity of such a pulsar is a direct manifestation of the quantum origin of the 
axionic siren. 
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We should mention that the above discussion ignores a number of important physical effects, 
such as the back reaction from self-interactions in the axion cloud, and one graviton annihilation 
of two axions in the same level, but they do not change the qualitative picture presented above. 
These effects will be discussed in an upcoming paper [71] . 

2.5.3 Superradiance and the QCD Axion 

For f a ~ M p i the QCD axion mass is 5 x 10~ 13 eV; given that the smallest stellar size black 
holes may have masses down to ~ 2M , their spin measurements may probe axion masses up to 
3 x 1CT 10 eV, exploring the parameter space of the QCD axion deep in the anthropic regime. In 
terms of the QCD axion decay constant this translates into 

f a > 2 x 10 16 GeV, (86) 

i.e. down to the GUT scale — the natural scale for the PQ symmetry breaking. Note that this 
reach is in a region where all electromagnetic couplings are suppressed and is also independent 
of the axion's cosmological abundance. From the already existing spinning black hole candidates, 
the best bound is provided by LMC X-l — a 1OM black hole with spin a/R g ~ 0.91 [66]. These 
values suggest a bound for f a at the level 

f a < 2 x 10 17 GeV , 

which is significantly below the Planck scale. 

Again, both these numbers should be considered as indicative (and conservative) estimates. 
Our preliminary results including the effects mentioned in the previous section and not accounted 
here, indicate that the actual bounds are likely to be stronger by a factor ~ 2. 

To conclude, it is worth noting that one could use the superradiance effect to place bounds on 
the photon and graviton masses. In the photon case, though, an electromagnetic field resulting 
from the instability will interact strongly with the matter in the black hole accretion disk and is 
most likely to dissipate. For the graviton, only the largest galactic black holes have a chance to 
compete with a limit from the binary pulsar timing. It was already suggested that these can be 
used to probe massive gravity models by gravitational wave measurements sensitive to the presence 
of black hole hair [72]; in order to exploit the superradiance effect it is necessary to observationally 
confirm the presence of an ergosphere region, which is not automatic in massive gravity, where the 
Kerr metric gets modified. An interesting consequence of a non-vanishing graviton mass would be 
that in the heavy mass regime graviton emission becomes impossible, and does not prevent the 
graviton cloud from becoming an order one perturbation. 

3 Discussion and Future Directions 

We see that taking seriously the QCD axion as a solution to the strong CP problem together with 
generic properties of axions in string theory leads to a rather unconventional set of predictions for 
forthcoming cosmological and astronomical observations. Of course, it would be exciting to observe 
any of the signatures discussed in the current paper — rotation of the CMB polarization, a step in 



32 



the CDM power spectrum or a gap in the mass spectrum of rotating black holes. However, a really 
distinctive property of the presented scenario is the expectation of many light axion fields. This 
multiplicity can reveal itself by giving rise to all of the above signals simultaneously. Even more 
interestingly, precision measurements of the CDM power spectrum and spectroscopy/gravity wave 
signals from rotating black holes are capable not only of providing the evidence for the presence 
of many light axions but also to measure values for many parameters of the same physical origin, 
such as axion masses/PQ breaking scales, and initial misalignments angles. With large enough 
number of axions this may allow to make quantitative tests of different scenarios for the statistical 
distribution of these parameters in the landscape and/or of inflationary measures. 

As discussed in the Introduction, if many ultra-light scalar fields were to be observed, it is 
tempting to suggest that they may be responsible for the scanning of the vacuum energy — at least 
at low scales. Indeed, in the current framework there are at least three natural sources for the 
scanning. First, the vacuum energy can be scanned at the high scale by fluxes, as in the Bousso- 
Polchinski scenario. Second, we expect a large number of real scalar moduli <f>i, their potential is 
generated as a result of SUSY breaking and depends on the SUSY breaking vev F and the string 
scale M s in the following way, 

Vfa) = F* usy f(fa/M s ) . 

This potential may provide a possibility to scan the vacuum energy below the SUSY breaking 
scale. Note, that unlike in the toy field theory landscape of [31], the maximum scale of scanning 
in this case is parametrically above the mass scale F/M s for these fields. Finally, the scanning can 
be done by axions. In this case the masses of the fields are also parametrically below the overall 
scale of the potential. The difference with the scanning by real moduli is that the axion masses 
are exponentially sensitive to the parameters of compactification, and likely to be distributed over 
many orders of magnitude rather than being concentrated around one particular scale, giving rise 
to the situation close to the Wilsonian scanning. 

The latter possibility raises a number of theoretical and phenomenological issues. First, as 
we already discussed, the scanning is only possible, provided individual axion potentials have 
non-degenerate minima. This implies that axions responsible for the scanning are associated with 
a gauge sector strongly coupled at low energies. Note, that by itself this does not guarantee 
the presence of multiple axion vacua. For instance, the axion potential in pure gluodynamics is 
quadratic in the large N limit [73], so that there is a single axion vacuum. 

Note that in this example the strongly coupled sector itself, even in the absence of the axion, 
possess iV non-degenerate metastable vacua. These vacua are useless for the scanning, however. 
First, not all of them are extrema of the axion potential. Moreover, even if in a more complicated 
setup, this kind of vacua become local minima of the axion potential, they still do not allow 
to scan the vacuum energy if the strong coupling scale A of the gauge sector is lower then the 
expansion rate of the Universe during inflation, H in f. Indeed, at the moment when at the FRW 
stage the expansion rate drops below A, the gauge sector experiences the (de) confinement phase 
transition (assuming it was not reheated, otherwise it happens later) and a network of domain walls 
separating different vacua is formed. Afterwards, bubbles of the lowest energy vacuum expand at 
the expense of the other vacua and the Universe always ends up being in the lowest energy state. 
This has to be contrasted with the axion vacua that exist only in the presence of a dynamical 
axion. Those are separated by a large distance of order f a in the field space. So if H in f <C f a 
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(but, possibly, H in f ^> A) the axion fluctuations are negligible compared to the distance between 
vacua and the network of domain walls does not form. We see that an important theoretical issue 
that has to be understood is what are the conditions for developing non-degenerate vacua in the 
axion potential and how generic such a situation is. 

If axions scan a significant fraction of the vacuum energy in the Wilsonian way, there should 
be several hidden gauge sectors per decade of energy at low energy scales. How is it possible that 
they avoided being detected so far? The natural answer to this question would be that these 
sectors are well separated from us along extra dimensions, so that locality sequesters them away. 
For instance, if inflation were due to some well-localized process in extra-dimensions (such as the 
brane inflation) this would explain why only the visible sector was reheated. 

On the other hand, for some purposes a large separation along extra dimensions may not be 
enough to realize sequestering at the desired level. For instance, it is conceivable that some of 
the hidden sectors posses massless U(l) gauge factors and have light (with masses smaller than 
~ 10 keV) fermions charged under these gauge bosons. In this case, there are extremely strong 
bounds from star cooling on the kinetic mixing coefficient e between the photon and an extra 
massless gauge boson e < 1CT 14 [H]. On the other hand, it was argued [75] that even if the 
extra U(l) is separated from our sector by a large distance in string units, this mixing can be 
mediated by light (compared to the string scale) closed string modes at the unacceptable level. 
Consequently, the smallness of such kinetic mixings requires a separate explanation. For instance, 
it might be a consequence of a large mass for the closed string modes capable of mediating such a 
mixing. It would be enough if such a mass were generated just in the vicinity of where the photon 
is localized. 

Note, that similar to the logic in the Introduction that led us to the expectation that many light 
axions may be present, one can also make the case for extra hidden gauge sector^] Consequently, 
these phenomenological issues remain even if axions are not responsible for the scanning. 

Another feature that can be generic in the compactification manifold is the existence of warped 
throats [76], W7\ 178] which generate hierarchically small mass scales by the Randall-Sundrum mech- 
anism [79]. Warping can naturally lower the value of f a below Mqut ~~ fa can now vary over a 
large range of mass scales, affecting the phenomenology of these axions. For example, the pre- 
dicted suppression in matter perturbations due to light axions changes, or, as explained in Section 



2.5.3[ the effect of superradiance can be inhibited due to the larger possible coupling to E ■ B. 
In addition, the scale of gauge sectors in warped throats can naturally be well below the string 
scale and may be giving rise to new ultra-light fields. These ultra-light fields can also include 
higher spin excitations gravitationally coupled to us. The phenomenology of axions with varying 
f a below the GUT scale, as well as possible implications from the existence of ultra light hidden 
sectors will be the subject of a companion paper. 

4 One difference is that unlike for the QCD axion it is possible to argue that the existence of at least one gauge 
sector is required anthropically. 
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5 Appendix: Superradiance as Anti-tunneling 

In Section |2.4| we have introduced the description of the exponential superradiant instability as 



a tunneling phenomenon, based on the Schrodinger problem (62), (63). Here we will discuss a 



simplified version of the problem which clearly demonstrates how the superradiance sets in. Let 
us for this purpose consider the ID Schrodinger problem 



tf" + (u 2 -V(x,uj\v = 



(87) 



where the potential is given by (see Figure [9]) 
V(x,uj) 



oo . 



a5(x + L) - 2mu + LU - m 2 oj\ 1 - 0(x + L) 



x > : 
x < 



(88) 



in the well of (88); and the 



m 2 u;i on the horizon's side of the 5-barrier, which is chosen to encode 



The potential (88) faithfully represents (63); it has: a mirror at large distances from the black 
holes, here modeled by the infinite potential barrier at the origin; the centrifugal potential at 
distances of the order of black hole's gravitational radius, represented by the strongly repulsive 
(^-function at x — —L, with a/u ^> 1; the potential well in between, with a negligible influence on 
the modes near the top of the well, reflected in our choice of V 
potential V = 2muo + uo 

the dispersion relation of a wave with angular momentum m in the near- horizon limit, where 
the wavevectors obey k 2 = (u — muj + ) 2 . Note, that albeit this problem is straightforward to 
solve, it is not of the usual self-adjoint form as commonly encountered in quantum mechanics on 
Hilbert space. Instead, the potential itself depends on the eigenvalue, which as we will now see 
has dramatic consequences for the existence of the instability. 

Let us now solve this equation. Away from the barriers, the solutions are given by linear 
combinations of free waves e ±iqx , with q = \Juj 2 — V in the corresponding region. At x = 0, 
where the wavefunction enters the barrier, we set \1/ = since the barrier is infinite both in height 
and in width, and penetration depth therefore vanishes. Far to the left of the ^-barrier, as we 



noted in Section 2.4, the regularity at the horizon picks the wave with the momentum dependence 
as the correct asymptotic form of the eigenmodes of (62), (63) [SB], where k 



exp^ 



-%k A 



x 



is 
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Figure 9: A toy model potential encapsulating all the salient features of the superradiant insta- 
bility: infinite potential barrier, representing the mass mirror; a potential well with V = 0, a 
reasonable approximation for the modes near its top; a strongly repulsive Dirac 5-function po- 
tential simulating the centrifugal barrier; and a potential V = 2muj + uj — m 2 ^, modeling the 
dispersive properties of the waves far in the ergoregion as they approach the horizon, where they 
satisfy k 2 = (u — muj + ) 2 , such that the net potential there is large and positive for large eigenval- 
ues, (u/ on the figure), but is negative for superradiant modes (denoted by uj on the figure). The 
key is that the superradiant modes have larger momentum outside the well than inside it. 



formally the positive root of q 2 = uj 2 — V = (uj — muj + ) 2 ; thus k + = uj — mu) + . Note, that when 
uj < muj + , this root is negative: it describes a wave which is incident on the 5-barrier from the 



left. Nevertheless, the group velocity of this wave is still positive, v g = 4r = 1, implying that any 
wavepacket composed of such modes moves away from the 5-barrier, even if the phase velocity of 
the wave changed sign. This is of course the first hint of the superradiant behavior in our simple 
toy model problem: the waves carrying away negative energy to the left leak into the potential well 
on the right, and amplify the state bound up in there. In fact, for the modes < uj < muj + we see 
that this boundary condition really represents the exact opposite of the usual textbook examples 
of tunneling, since the 'free' wave impinges into the barrier, and so we can dub it 'anti-tunneling'. 



What remains is to complete the determination of the spectrum of (87), (88) by matching the 



waves in the well and out in the ergoregion across the 5-barrier. There, the solutions obey 

^ L \ x= - L = Vr\ x= - L , ^'r\x=-l - ^'l\x=-l = a*U=-L • (89) 
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Substituting 

L = Be- i{uJ - mw+)x , ^r = A sm(ux) , (90) 

where these functions are picked to satisfy the boundary conditions far to the left and at the 
origin, respectively, we get a secular equation for the eigenvalues u: 

cu cot(u>L) + a = i(uj — mu+^j . (91) 



Two comments are in order here. First: the secular equation (91) should come as no surprise, as 



we are dealing with a problem that involves three modes and four boundary conditions, thanks 



to the horizon ingoing condition and infinite barrier cutoff enforcing (90), and (89) relating their 
parameters at the 5-barrier. Physically, this merely means that the potential well of Figure [9] 
can only accommodate states which have right phases to fit inside it, as it enforces two-sided 
boundary constraints. This is expected, as any bound state spectrum is indeed discrete. Second: 
the eigenfrequencies all have imaginary parts. This is because the 5-barrier is not impenetrable, 
unlike the infinite barrier at the origin, but allows leakage which links the bound states with 
the continuum of outgoing waves to the left. This is indeed familiar from the usual tunneling 
problems, and reflects the fact that under time evolution the bound states can evolve into free 
states that escape to infinity. In the standard tunneling problems, however, the imaginary parts 
of the eigenfrequencies are all negative, which means that the bound states decay as time goes on. 
That is again a natural consequence of a setup which one adopts, which is that a quantum state 
is prepared in a surrounding vacuum and allowed to fall apart by the emission of a wave through 
a partially transparent barrier. 



The unusual feature of (91) is that many of the eigenmodes have positive imaginary parts. 
This can be readily seen by setting 

to — u r — iT , (92) 



and solving (91 ) in the limit luL > 1 (which means that the well has many closely populated bound 
states, but that we look for ones near the top) and a/u ^> 1 (so that the 5-barrier is strongly 
repulsive, and the bound states are long-lived, allowing for an adiabatic approach to the problem). 
These approximations directly follow from the black hole instability considerations, which we are 



principally interested in. Then, the secular equation (91) splits into two 



u r cot(uj r L) + a + 0(T) =0, (93) 
T cot(w r L) = -(uj r - muo + ) + 0(Tuo r L) . (94) 

When a/u ^> 1 the standard magic of cotangents in the tunneling eigenequations comes to 



the rescue: the eigenmode real parts are determined by the limiting form of the equation (93), 
cot(avL) — —a/u r , and are approximately 

mr , , 

u, r ~— + .... (95) 



This implies that the solutions of the other secular equation (94) are 



a 
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This is the key equation for understanding the superradiant instability. We see immediately 



that the modes in the regime < u r < muj + have T < 0, which from Eq. (92) implies that 
these modes have positive imaginary contributions to the eigenfrequency. This means that, as 
time goes on, these modes grow exponentially. They are unstable, absorbing energy and spin from 
the black hole, rather than the other way around. For them, indeed, the 'transmitted' wave is 
really incident as k + = u — mw + < 0, and so they anti-tunnel, as we anticipated above. Finally, 
we also see that the instability scale T peaks at uo T ~ mw + /2, as that is the extremal value of V 



in Eq. (96). This of course presumes that the 5-barrier is independent of the eigenfrequency cu, 
as well as that the mirror barrier is fully impenetrable, as modeled by our infinite barrier at the 
origin. In the real black hole problem, neither is true. The real centrifugal barrier depends on 
the eigenmode, through the respective value of I. And the mirror is not completely impenetrable, 
being eventually completely transparent to the modes with u > fi. All of this will modify the 
precise formulas for the superradiant eigenfrequencies and the instability time scales. Nonetheless, 
we believe the above toy model is instructive and accurately reproduces the essential features of 
the instability. 
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